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ABSTRACT 


This  report  documents  the  procedure  for  generating  a  "real  world" 
ephemeris  tape  to  be  used  by  The  Aerospace  Corporation  and  selected  con¬ 
tractors  in  Phase  0  of  the  Autonomous  Navigation  System  (ANS)  contract. 

A  "real  world"  geopotential  model  was  developed  by  modifying  a  state-of- 
the-art  reference  geopotential,  using  Kaula's  degree  variances  as  a  guide. 

To  obtain  "real  world"  atmosphere  daia,  the  acceleration  profile  experienced 
by  the  first  in  a  recent  series  of  low  altitude  satellites  to  have  an  on-board 
low-g  accelerometer  was  suitably  scaled.  The  resulting  ephemeris  is 
displayed  and  subjected  to  various  consistency  tests. 
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SECTION  I 


INTRODUCTION 


A  "real  world"  ephemeris  tape  has  been  generated  for  use  by  The 
Aerospace  Corporation  and  selected  contractors  in  Phase  0  of  the  Autonomous 
Navigation  System  (ANS)  contract.  This  ephemeris  is  based  on  an  analytic 
geopotential  force  model  and  an  accelerometer  data  tape  (in  lieu  of  an  atmo¬ 
sphere  model) . 

To  serve  as  a  reference  and  a  standard  for  comparison,  a  "model” 
ephemeris  was  also  generated.  This  model  ephemeris  employed  the  same 
initial  .onditions  that  were  to  be  used  in  the  "real  world”  ephemeris  and  used 
force  models  that  are  representative  of  the  state  of  the  art. 

Independent  of  the  model  ephemeris  generation,  a  geopotcnliai  model 
which  simulates  the  real  world  was  developed.  To  represent  the  atmosphere, 
a  span  of  accelerometer  data  was  taken  from  an  actual  satellite  flight  whose 
orbit  was  quite  similar  to  that  chosen  for  the  reference.  These  accelerations 
were  scaled  to  meet  certain  criteria  to  make  them  suitable  for  use  as  a  real 
world  atmosphere  representation.  The  detailed  development  of  these  models 
is  described  in  the  billowing  sections. 

Comparisons  between  the  model  and  real  world  ephemeridos  were  per¬ 
formed  to  ensure  that  the  real  world  case  satisfied  certain  consistency  tests. 

Figure  ^-1  shows  diagrammaticaliy  the  tasks  described  above.  The 
work  on  each  task  is  documented  in  the  section  noted  parenthetically  within 
the  box  describing  the  task. 
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SECTION  II 


INITIAL  CONDITIONS 


The  low  altitude  reference  orbit  for  the  ANS  has  initial  perigee  occurring 
within  on^-half  hour  of  local  noon  at  approximately  45°N  latitude.  The  initiaL 
apogee  altitude  was  200  n  mi  and  perigee  altitude  was  80  n  ini  above  a  spherical 
earth;  its  period  was  about  89.7  min. 

Epoch  time  for  this  orbit  was  chosen  as  2230  GMT  on  29  February  1980. 
The  initial  conditions  are  shown  in  Table  II- 1  in  four  coordinate  representations. 
These  are  conditions  at  an  ascending  node  for  an  orbit  having  the  previously 
described  charact  eristics. 

The  first  set  is  an  earth-centered  inertial  (ECI)  set  which  has  as  its 
reference  plane  the  true  equator  at  epoch  and  as  its  reference  direction  the 
mean  equinox  at  0  hr  GMT  of  the  date  of  epoch. 

The  second  set  is  an  inertial  spherical  coordinate  set  having  o  and  6  as 
the  right  ascension  (measured  from  the  X  axis,  positive  eastward)  and  geo¬ 
centric  latitude  (declination)  ol  the  vehicle;  p  is  the  angle  between  the  velocitv 
vector  and  the  geocentric  vertical  (the  vertical  flight  path  angle);  Ac  is  *he 
azimuth  of  the  velocity  vector  from  true  north,  measured  eastward  on  a  plane 
normal  to  the  geocentric  vertical;  ar.d  R  and  V  are  the  respective  magnitudes 
of  the  position  and  velocity  vectors. 

The  third  set  is  an  earth-fixed  spherical  coordinate  set  having  X.  and  6 
as  the  geodetic  latitude  and  longitude  of  the  vehicle;  y  as  the  angle  between  the 
velocity  vector  and  the  geocentric  horizontal  (the  horizontal  flight  path  angle); 
and  Az,  R,  and  V  as  nreviously  defined. 

The  fourth  set  is  the  classical  element  set,  where  a  and  e  arc  the  semi¬ 
major  axis  and  eccentricty,  respectively;  i  is  the  inclination  of  the  orbit  plane; 

S7  is  the  right  ascension  of  the  ascendirg  node;  w  is  the  argument  of  perigee,  the 
angle  between  the  direction  of  perigee,  and  the  line  of  nodes  measured  from 
ascending  node  to  perigee;  and  r  is  a  reference  time  indicating  the  time  ol  last 
perigee  passage  in  minutes  from  epoch. 
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Table  II- 1 .  Four  Representations  of  Initial  Conditions 


(  1)  Cartesian  (ECI) 

(2)  Spherical 

X 

-  1.  59387494563E  L7  ft 

a 

136.  390753744  deg 

Y 

1. 518517  14534E+7  ft 

6 

1 . 0E-20  deg 

/. 

3. 8419943  1557E-  15  ft 

p 

90. 7298014230  deg 

X 

6.  1 6441  28293 5E  +3  fps 

Az 

340  deg 

Y 

6.  00674668844E  +3  fps 

R 

2.  2013005 9  195E +7  ft 

/. 

2.  363  1  2409387E  +4  fps 

V 

2.  514<i8815029E  M  fps 

(3)  Ea  rth- Fixed  Sphe  rical 

(4)  Classical  Elements 

\ 

0  d  e  g 

a 

2.  177 63 946Q  1  6E  +7  ft 

o 

360  deg 

e 

1.  67413 647743 E- 2 

7 

-0.  7 2980  14  23  0  deg 

i 

1  1 0  deg 

A/. 

340  deg 

n 

130.390753744  deg 

i\ 

2.  20  1  3 00  59  1 95E  +7  ft 

u; 

13  1.  14363965  1  deg 

V 

2.  5  14988  15029E+4  fps 

T 

-57.36926627  min 

The  Greenwich  hour  u.ngle  at  0  hr  GMT  of  the  assumed  epoch  data  is 
157.9667413  deg.  This  value  is  needed  to  facilitate  transformations  among 
these  four  representations. 

The  reference  ballistic  coefficient  (C^A/W)  for  this  vehicle  is  0.  02  ft 
This  value  is  considered  typical  oi  future  high  drag  satellite  system 
configurations. 


/lb. 
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SECTION  HI 


"MODEL"  EPHEMERIS 


The  TRACE  computer  program  (Ref.  1  and  2)  was  utilized  in  generating 
the  model  ephemeris.  The  initial  conditions  previously  discussed  were  used 
to  generate  an  ephemeris  for  20  orbits  spanning  approximately  1-1/4  days. 

Force  models  used  in  the  model  ephemeris  included  a  6th  degree  and 
order  versicn  of  a  Guier  8th  degree  and  order  geopotential  (Ref.  3),  and  the 
Jacchia-Walker-Bruce  atmosphere  model  (Ref.  4).  A  90-day  mean  of  the 
10.7  cm  solar  flux  F^q,  denoted  Fjq,  is  used  as  a  reference  about  which  the 
short-term  decametric  flux  is  measured.  Values  of  F  q  =  F^q  =  220  w/M^/cps 
were  used  in  vhe  model  ephemeris.  A  value  of  the  planetary  magnetic  index 
ap  =  20  was  also  used.  Appendix  A  contains  a  series  of  plots  showing  the 
ground  trace  of  this  veh'cle  as  a  function  of  time.  For  reference  and  for  later 
comparison  purposes,  Table  III  - 1  gives  the  times  and  ECI  coordinates  of  the 
vehicle  at  each  ascending  node. 


Preceding  page  blank 
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Table  Ill- 1 .  Nodal  Conditions  for  "Model"  Ephemeris 


Rev 

No . 

Time  of  Ascension 
Nodal  Cross 
/Mo/  Day/ Y  r  \ 

VHr /Min/ Sec/ 

Position  (ECI)  XYZ 
(ft) 

■ . —  - . . 

Velocity  (ECI)  XYZ 
(fps) 

0 

2/29/80 

22/30/0.  00 

-  1.  59387495E+7 

1.  5183  1715E+7 

3. 84 19943  2E-  15 

6.  16441283E+3 

6.  00674669E+3' 

2.  363  12409E+4 

< 

l 

2/29/80 

23/59/34. 88980 

-  1. 59849728E+7 

1.  51290676E+7 
-3. 42450058E-4 

6.  145263  10E+3 

6.  02746746E+3 

2.  36342563E+4 

2 

3/1/80 

1/29/9. 04968 

-  1. 60306694E+7 

1.  5 07 4 9 26 IE +7 
-1.  27921077E-3 

6.  1  258253  2E+3 

6.  0481 877 2E  +3 

2.  36374575E+4 

3 

3/1/80 

2/fc  8/42.  38501 

-  1 . 60763  27  9E +7 

1.  50203  121E+7 
-3. 933447 1 8E-4 

6.  105953  10E+3 

6. 06958938E+3 

2.  36408159E+3 

4 

3/1/80 

4/16/22.  60609 

-1. 612223 26E+7 

1 . 4965  1661E+7 
-3. 2 1 860488E-4 

6.  08579606E+3 

6.  09047463E+3 

2.  36442966E+4 

5 

3/1/80 

5/57/46. 87890 

-  1 . 61 677  347  E  f7 

1 . 49098307 E  +7 

-  1.  2969 1 562E-3 

6.  066 1 6627E+3 

6.  1  1094373E+3 

2.  3  647 844 1 E +4 

6 

3/1/80 

7/27/18.  02516 

-  1. 62126068E+7 

1 . 4854606 1E+7 
-5. 61057958E-5 

6,  04703706E+3 

6.  13214743E  (-3 

2.  3652798E+4 

7 

3/1/80 

8/56/48.31871 

-1. 625744 1  2E  +7 

1. 47994091E+7 
-  1 , 34453  140E-3 

C.  027  63290E+3 

6.  15356280E+3 

2.  36544996E+4 

8 

3/1/80 

10/26/ 17. 81324 

-1.63022491E+7 
1.47441290E+7 
-4. 843 07  36 2E- 5 

6.  008147  12E+3 

6.  1747  1305E+3 

2.  3 657 5 294 E +4 

9 

3/1/80 

11/55/46.47345 

-  1. 63468743E+7 

1 . 468896 14E  +7 

-  1. 26815078E-3 

5. 98856274E+3 

6.  19570106E+3 

2.  36603321E+4 
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Table  III—  1 .  Nodal  Conditions  for  "Model"  Ephemeris  (Continued) 


Rev 

No. 

Time  of  Ascension 
Nodal  Cross 
/Mo/  Day  /  Yr\ 
\Hr/Min/ Sec/ 

Position  (ECI)  XYZ 
(ft) 

Velocity  (ECI)  XYZ 
(fps) 

10 

3/1/80 

13/25/14.30215 

-  1.63916332E+7 
1.46338136E+7 
-4.  1 800963  9E-4 

5.  968337 27E+3 

6.  2 1 60867  3E+3 

2.  36630337E+4 

11 

3/1/80 

14/54/41.38993 

-  1.  64363722E+7 

1 . 457  83  174E+7 
-2.  46736766E-4 

5.  947 5 2069E+3 

6.  23585677E+3 

2.  36659569E+4 

12 

3/1/80 

16/24/7. 75793 

-  1. 64807430E+7 

1 . 45224  105E+7 

-  1. 30417512E-3 

5.  92644969E+3 

6.  25589753E+3 

2.  36690683E+4 

13 

3/1/80 

17/53/33.34407 

-  1.  65248655E+7 
1.44662676E+7 
-6.31 984487E-4 

5.  90592687E+3 

6.  27 676825E  +3 

2.  36720090E+4 

14 

3/1/80 

12/22/58. 06510 

-  1 . 65687 876E  +7 

1.441 0066 1 E+7 
-6.  78497093E-5 

5.  88630 149E  33 

6.  29783075E4-3 

2.  3 6747 293E+4 

15 

3/1/80 

20/52/21. 87667 

-  1. 66124592E+7 

1. 43537098E+7 
-3.  2273460 IE-4 

5.  86662222E  33 

6.  3  1863 64 2E +3 

2.  3677  5  1 53  E  34 

16 

3/1/80 

22/21/44. 78377 

-  1. 66559668E+7 

1. 42970491E+7 
-1. 05750460E-3 

5.  84 64 24 5 9E +3 

6.  33Q35065E+3 

2.  36805229E+4 

17 

3/1/80 

23/51/6. 84696 

-  1. 66991  600E+7 

1. 42403605E+7 

-  1. 301527 80E-3 

5.  825^777  1E+3 

6.  3  6O45009E  1-3 

2.  36837 258E+4 

18 

3/2/80 

1/20/28.  09812 

-  1. 67418053E+7 
1.41837  1 85E  +7 

-  1.  32891  135E-3 

5.  80526738E  33 

6.  37930109E+3 

2.  3  687  1  191E+4 

19 

3/2/80 

2/49/48.44448 

-  1. 67842644E+7 

1.41  2667  81  E+7 

-  1.  32615034E-3 

5. 784  1 1  142E+3 

6.  39985«57E+3 

2.  3690683 5E  34 

20 

3/2/80 

4/19/7.90466 

-  1. 68269602E+7 

1 . 40690575E+7 

-  1.  30997504E-3 

5.76259368E+3 

6.  4201  2503E+3 

2.  ?  694  37  83  E  34 
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SECTION  IV 


"REAL  WORLD"  GEOPOTENTIAL 


A.  SIMULATED  GRAVITY  FIELD  GENERATION 

For  the  real  world  case,  the  geopotential  was  expressed  as  a  spherical 
harmonic  expansion  of  the  following  form: 


U 


GM 


-t(W 

x\~l. 


n 


Cno?n(sin  +  ]CPnm(sin  *)(C 


cos  m\  +  S  sin  ml) 
nm  nm 


m=l 


(IV- 1) 


where 

U  =  geopotential 

GM  =  gravitational  constant  times  the  mass  of  the  earth 
=  equatorial  radius  of  the  earth 
r  =  radial  distance  from  the  earth's  center  of  mass 


\ 


pm 

n 

C  ,  S 
nm  nm 


geocentric  latitude 
geocentric  longitude 

normalized  associated  Legendre  functions 

normalized  cosine  and  sine  coefficients  of  the  nth  degree 
and  mth  order 


When  the  gravity  field  is  expressed  as  such  a  spherical  harmonic  expansion, 
it  is  completely  described  by  the  choice  of  spherical  harmonic  coefficients 
(C,  S).  This  study  thus  required  that  a  suitable  set  of  C's  and  S's  be  selected. 

The  set  of  spherical  harmonic  coefficients  used  for  this  study  was 
generated  as  follows.  First,  a  set  of  C's  and  S's  was  selected  from  one  of  the 
currently  available  geopotentials  published  in  the  literature.  Second, 


Preceding  page  blank 

1 1 


Table  IV- 1.  Gravity  Degree  Variances 


Degree 

Mean  of  Coefficients 
with  Standard  Deviation 
About  the  Mean 
of  Reference 

Mean  of  Perturbations 
with  Standard  Deviation 
About  Mean  X  106 

2 

-96.6  ±  216.6 

0. 0015  ±  0. 036 

3 

0.  798  ±  0.  853 

0. 0072  ±  0. 019 

4 

0.  154  ±  0.  525 

0.  0123  ±  0.  120 

5 

-0.094  ±  0.361 

-0.028  ±  0. 075 

6 

-0. 152  ±  0. 216 

-0. 023  ±  0. 054 

7 

0.  026  ±  0.  192 

0.0018  ±  0. 045 

8 

0.  026  ±  0.  13  1 

-0. 0069  ±  0. 036 

9 

0.013  ±  0.  118 

0. 0055  ±  0. 042 

10 

-0. 023  ±  0.095 

-0. 009  ±  0. 035 

1 1 

-0. 036  ±  0. 076 

0.  014  ±  0.  027 

1Z 

-0. 0014  ±  0.  109 

0.  0039  ±  0.  033 

13 

0. 015  ±  0. 07 V 

0.  0083  ±  0.  028 

14 

0.  015  ±  0.  081 

-0. 0036  ±  0.  029 

.15 

0.  0081  ±  0.  063 

0. 003  ±  0.  032 

lb 

-0. 0151  ±  0. 0685 

-0.  0095  ±  0.  0278 

17 

-0.  0032  ±  0.  0539 

-0. 0003  ±  0.  0268 

18 

-0. 0043  ±  0. 060 

0.  0006  ±  0. 030 

19 

0.  0036  ±  0.  051 

-0. 0022  ±  0. 029 

ZO 

-0. 00066  ±  0. 025 

0.0037  ±  0. 022 

Z1 

-0. 0020  ±  0. 01 1 

0. 00052  ±  0. 02b 

Z2 

-0. 0011  ±  0. 015 

0. 00012  ±  0. 023 

23 

0. 00055  ±  0. 023 

0. 00284  ±  0. 020 

24 

0.  0027  ±  0.  010 

-0. 0012  ±  0. 021 

perturbations  were  then  algebraically  added  to  the  set  of  coefficients  in  a  way 
that  left  the  lower  degree  features  of  the  field  unchanged  but  introduced 
increasing  amounts  of  variation  with  increasing  degree.  Third,  additional 
coefficients  were  generated  for  those  terms  not  provided  in  the  original  set 
to  increase  the  degree  of  the  resulting  simulated  gravity  field  to  the  maximum 
size  permitted  by  existing  software. 

Perturbations  to  the  reference  field  were  generated  in  the  following 
manner:  First,  a  degree  variance  was  chosen  for  every  degree  2  through  24. 
Table  IV- 1  shows  the  degree  variance  used  in  generating  this  field,  the  values 
obtained  and  published  by  W.  M.  Kaula  (Ref.  5),  and  the  values  ublained  from 
the  reference  set.  As  can  be  seen  from  the  table,  the  degree  variances  were 
chosen  to  introduce  very  small  perturbations  in  the  lower  degree  field  with  an 
increasing  variation  as  upper  degree  fields  were  generated,  until  finally  at 
19th  degree  a  field  would  be  produced  that  roughly  obeyed  the  properties 
reported  by  Kaula. 

For  every  degree  variance  of  Table  IV - 1  (column  1),  a  set  of  coefficient 

perturbations,  AC  ,  AS  ,  was  chosen  from  a  random  number  generator 
nm  nm 

such  that  the  sum  of  the  squares  of  these  coefficients  would  satisfy  the  equation 


2 

a 


n 


1 


n 

E 

m=0 


_ 2  _ 2 

AC  +  AS 
nm  nm 


(IV-2) 


where  is  the  degree  variance  and  g^  is  the  mean  equatorial  gravity. 

Finally,  the  coefficients  generated  by  this  process  were  algebraically 
added  to  the  coefficients  of  identical  degree  and  order  in  the  reference  set, 
and  the  result  was  then  used  as  the  final  representation  of  the  earth's  gravity 
field. 

Using  the  above  method,  a  field  was  generated  complete  to  24th  degree 
and  order.  This  field  above  19th  degree  is,  essentially,  a  random  field 
conforming  only  to  the  constraint  of  Eq.  (IV-2).  For  19th  degree  and  under, 
the  field  is  a  slightly  varied  form  of  the  reference  field.  This  can  be  seen 


13 


more  clearly  in  Table  IV-2,  which  gives  the  results  of  averaging  both  the 
reference  set  and  the  perturbation  set  of  coefficients  over  all  coefficients  of 
the  same  degree.  The  sigma  presented  in  Table  IV -2  is  the  variance  of  the 
coefficients  of  a  particular  degree  about  the  mean  of  the  coefficients  of  that 
same  degree.  This  sigma  should  not  be  interpreted  as  any  measure  of 
accuracy  of  the  coefficients  but  is  included  to  provide  some  measure  of  the 
variation  of  coefficients  of  a  particular  degree. 

Table  IV-3  presents  the  final  set  of  coefficients  generated  for  this  study. 
These  coefficients  should  be  used  with  the  expression  for  earth  gravity  potential 
given  in  Eq.  (IV -1). 

B.  TEST  OF  SIMULATED  GRAVITY  FIELD 

The  field  specified  by  Table  IV-3  was  tested  for  its  effect  upon  the 
nominal  orbit.  This  test  was  conducted  by  performing  trajectory  difference 
runs  with  TRACE  (see  Ref.  1).  These  difference  runs  were  performed  by 
taking  the  nominal  epoch  vector  from  Table  II- 1  as  an  initial  condition  from 
which  an  orbit  was  obtained  by  numerical  integration  of  the  gravity  force 
model.  The  same  process  was  then  repeated  with  the  perturbed  gravity  field. 
The  difference  between  the  two  orbits  at  identical  times  was  then  resolved  into 
in-track,radial,  and  cross -track  differences. 

First,  a  difference  run  was  made  between  the  Guier  n  =  6,  m  =  6  model 
("model"  geopotential)  and  the  reference  model.  No  atmosphere  model  was  in¬ 
cluded.  The  purpose  of  this  run  was  to  establish  a  baseline  for  the  behavior  of 

a  high-degree  gravity  field.  (The  Guier  6.  6  geopotential  uses  a  value  of 

3  2 

GM  =  1. 4076538841E+  16  ft  /sec  ,  while  the  reference  model  has  a 
GM  =  1. 4076468597E+16  ft^/sec^.  )  Figures  IV-la,  b,  and  c  show  the 
behavior  of  the  in-track,  radial,  and  cross -track  differences.  The  in-track 
difference  shows  a  secular  growth  component  of  approximately  770  ft/rev 
together  with  a  periodic  component  with  a  maximum  excursion  of  1500  ft;  the 
radial  difference  shows  a  periodic  variation  with  a  maximum  amplitude  of 
approximately  700  ft;  the  cross-track  difference  has  a  maximum  periodic 
difference  of  1300  ft. 
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Table  IV-3.  "Real  World"  Normalized  Geopotential  Spherical  Harmonic  Coefficients  (Continued) 
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Table  IV-3.  "Real  World"  Normalized  Ceopotential  Spherical  Harmonic  Cojfficien;s  (Continued) 
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Table  IV-3.  "Real  World"  Normalized  Geopotential  Spherical  Harmonic  Coefficients  (Continued) 
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Next  a  difference  run  was  made  between  the  reference  model  and  the 
model  defined  by  Table  IV-3,  the  simulated  gravity  field.  These  in-track, 
radial,  and  cross -track  differences  are  presented  in  Figures  IV-2a,  b,  and  c. 
The  in-track  difference  shows  a  secular  growth  of  approximately  280  ft/rev 
plus  two  periodic  components:  one  with  a  period  of  1  rev  and  an  amplitude 
ri  approximately  800  ft  and  the  other  with  a  period  of  about  12  revs  and  an 
amplitude  of  1300  ft.  The  radial  difference  plot  shows  a  1  rev  periodic 
component  with  maximum  amplitude  of  approximately  7  00  ft,  as  does  the  cross¬ 
track  difference  plot.  As  before,  these  data  do  not  include  an  atmosphere 
force  model. 

It  was  concluded  from  these  results  that  the  simulated  gravity  field 
produces  orbit  perturbations  which  differ  significantly  from  those  of  the 
nominal  gravity  field  but  which  are  not  unrealistically  large.  The  results 
of  these  difference  runs  should  not  be  interpreted  as  estimates  of  ephemeris 
uncertainties,  since  by  suitable  orbit  determination  techniques  these  differences 
would  be  drastically  reduced  in  actual  practice.  The  plots  illustrate  one  of  the 
possible  shortcomings  of  this  simulated  model,  i.  e.  ,  the  field  is  too  smooth. 
Although  this  field  is  as  extensive  as  the  present  state  of  the  art,  there  is 
reason  to  believe  that  a  24th  degree  and  order  model  is  still  significantly 
smoother  than  that  experienced  by  a  vehicle  moving  in  the  earth's  gravity 
field. 
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Figure  IV-ia.  Intrack  Difference  Between  "Model"  ( 
and  Reference  Geopotential  (19,  19)  A 
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Figure  IV-ia.  Intrack  Difference  Between  "Model"  (6,  6)  Geopotential 
and  Reference  Geopotential  (19,  19)  Atmosphere 
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Figure  IV-ic.  Crosstrack  Diffe 
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Figure  IV-2b.  Radial  Difference  Between  Refer 
and  "Real  World"  Geopotentials 


Figure  IV-2b.  Radial  Difference  Between  Reference 
and  "Real  World"  Geopotentials 
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Figure  IV-2c.  Crosstrack  Difference  Between  R« 
and  "Real  World"  Geopotentials 


SECTION  V 


"REAL  WORLD"  DRAG  MODEL 


To  obtain  real  world  atmosphere  data,  the  acceleration  profile 
experienced  by  the  first  in  a  recent  series  of  low  altitude  satellites  to  have  an 
on-board  low-g  accelerometer  was  utilized.  The  orbital  elements  of  this 
density  measurement  system  (DMS)  vehicle  were  quite  close  to  those  of  the 
reference  orbit.  This  suggested  that  the  accelerometer  data  from  the  DMS  - 
suitably  scaled  -  could  be  used  as  an  input  real  world  drag  profile  to  the  ANS 
reference  vehicle.  Approximately  20  revs  of  continuous  DMS  data  were 
available  for  this  purpose. 

Figure  V-l  shows  a  comparison  of  the  model  drag  accelerations  and  the 
raw  DMS  sensed  accelerations.  It  is  readily  seen  that  their  amplitudes  do 
not  agree,  nor  do  points  of  maximum  d»*ag  in  the  two  sets  of  acceleration 
profiles  occur  at  the  same  time.  Approximately  240  sec  of  data  was  dubbed  on 
the  front  of  the  DMS  data  and  a  multiplier,  determined  by  the  ratio  of  the  total 
time  spans,  was  then  applied  to  the  result  to  force  alignment  of  the  points  of 
maximum  drag.  Since  the  orbital  periods  were  very  nearly  equal,  the 
multiplier  was  quite  close  to  1,  specifically  0.99770736.  This  resulted  in  the 
situation  depicted  in  Figure  V-2. 

If  this  data  were  used  directly  for  real  world  drag  acceleration  without 
further  scaling,  it  is  possible  that  the  resulting  integrated  real  world 
ephemeris  would  have  times  of  perigee  passage  out  of  phase  with  the  times  of 
maximum  drag  force.  Also,  the  magnitude  of  the  perigee  drag  acceleration 
itself  might  have  appeared  unreasonable  in  the  absence  of  any  scaling.  To 
prevent  this  from  occurring,  the  amplitude  of  the  data  was  scaled  so  the 
real  world  energy  loss  would  be  the  same  as  that  of  the  model.  Specifically, 
the  energy  loss  is  directly  proportional  to  the  velocity  loss 

AV  =  J" aD  dt 

where  a^  is  the  drag  acceleration. 
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Figure  V-2.  DMS  Dl 

Jl 

The  accelerometer  data  scaling  was  done  two  different  ways.  In  the 
generation  of  the  real  world  atmosphere  data  reported  on  herein,  the  energy 
loss,  or  AV,  was  preserved  on  a  rev-to-rev  basis.  A  second  candidate 
real  world  ephemeris  was  also  generated  in  which  the  total  20-rev  AV  was 
preserved  instead. 

The  slight  differences  between  the  two  real  world  ephemerides  using 

rev-to-rev  AV  scaling  and  20-rev  AV  scaling  were  considered  negligible  and 

the  former  case  was  chosen,  albeit  somewhat  arbitrarily,  to  be  the  real 

world  drag  profile  for  this  study.  Although  the  DMS  vehicle's  perigee  was 

2 

about  73  n  mi,  it  had  a  lower  ballistic  coefficient  (0.  008  ft  /lb)  than  that 
chosen  for  the  reference  vehicle  (0.  02  ft  /lb).  As  a  result,  the  rev-by-rev 
velocity  loss  for  the  DMS  flight  was  less  than  that  for  the  reference  orbit. 
Each  rev  of  the  DMS  data  was  scaled  by  the  ratio 

AV  loss  reference 
AV  loss  DMS 

which  is  a  number  greater  than  1.  Figure  V-3  compares  the  now-normalized 
DMS  data  with  the  model  drag  profile.  Each  rev  of  both  orbits  imparts 
‘'he  same  energy  loss  to  the  vehicle,  but  the  real  world  drag  at  any  given  point 
within  a  rev  generally  differs  from  that  of  the  model.  As  Figure  V-3  shows, 
it  would  be  extremely  difficult  to  reproduce  the  normalized  DMS  acceleration 
profile  with  an  atmosphere  model. 

The  normalized  DMS  drag  profile  was  written  on  tape  in  a  format 
suitable  for  use  as  an  accelerometer  input  to  TRACE. 
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Figure  V-3.  DMS  Data;  Scaled  Rev-by-Rev 


SECTION  VI 


"REAL  WORLD"  EPHEMERIS 


The  geopotential  described  in  Section  IV  and  the  accelerometer  drag 
profile  described  in  Section  V  were  used  as  force  models  together  with  the 
initial  conditions  described  in  Section  II.  The  gravitational  attractions  of  the 
sun  and  moon  were  also  included  to  make  the  real  world  ephemeris  as  near  to 
the  real  world  as  possible.  The  integration  was  carried  out  over  the  same 
20  revs  as  the  model  case. 

Since  the  real  world  drag  impulse  was  forced  to  agree  on  a  rev-to-rev 
basis  with  the  model  drag  impulse,  the  nodal  crossing  times  and  positions  for 
the  real  world  ephemeris  should  be  close  to  those  of  the  model  ephemeris. 
Table  Vl-1  contains  (in  the  same  format  as  Table  II-2)  these  conditions  for  the 
real  world  ephemeris. 

A  computer  listing  of  the  ephemeris  for  the  first  two  revs  is  included  in 
Appendix  B. 
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Table  VI- 1.  Nodal  Conditions  for  "Real  World"  Ephemeris 


Node 

Nodal  Crossing  Time 
/Mo/Day/Yr  \ 
\Hr/Min/Sec/ 

Position  XYZ 
(ft) 

Velocity  XYZ 
(fps) 

0 

2/29/80 

22/30/0. 0 

-1.  59387495E+7 

1.  5  1  83  17  1  5E+7 

3. 84 1 9943  2E-  1 5 

6.  16441283E+3 

6.  00674669E+3 

2.  36312409E+4 

1 

2/29/80 

23/59/34. 86323 

-  1.  59847865E+7 

1.  51293060E+7 
-4. 34333  50  IE-4 

6.  14523022E+3 

6.  02755269E+3 

2.  3 634 1 64 5 E +4 

2 

3/1/80 

1/29/9.  .07663 

-1. 60300679E+7 

1.  50750883E+7 
-9. 00100265E-4 

6.  1  2585377E+3 

6.  047 955  15E+3 

2.  3637 87 08E+4 

3 

3/1/80 

2/58/42.47044 

-  1. 60754250E+7 

1.  50202680E+7 
-3.7  257 1 3  21E-4 

6.  1 0597427  E +3 

6.  06962440E+3 

2.  364  1 5  889E+4 

4 

3/1/80 

4/28/15.08251 

-  1 . 6 1 214667 E +7 
1.496507  87E+7 

7.  1 3  800964E-  5 

6.  085647 28E  +3 

6.  09052957E+3 

2.  3  644  96 6  2E +4 

5 

3/  1/80 

5/57/47. 03614 

-  1. 61672548E+7 
1.49094396E+7 
- 1.45264228E-3 

6.  06617  102E+3 

6.  1  1093468E+3 

2.  36485783E+4 

6 

3/1/80 

7/27/18. 19338 

-  1. 621 22626E  +7 

1 .48540354E+7 

2.  97  1  127  82E-5 

6.  0464  8130E+3 

6.  1 3  2 1  87  03  E  +3 

2.  3652257  1  E+4 

7 

3/1/80 

8/56/48. 50253 

-  1. 6257  1392E+7 

1 . 47  988953  E  +7 

-  1.44689247E-3 

6.  027204  1  1E+3 

6.  15384  146E+3 

2.  3655  1«01  E+4 

8 

3/1/80 

10/26/ 17. 98143 

-  1 . 630 17 93 0E  +7 

1 . 4  7  4  3  4  6  2  5  E  +7 
-9. 23 92807 9 E- 5 

6.  007  54  997  E  +3 

6.  17  5 1947 3 E  +3 

2.  3  65  844  12E +4 

9 

3/1/80 

1 1/55/46. 64434 

-  1 . 63463 1 89 E  +7 

1 . 4688347  2E  L1 

-  1 . 6584  23 3  5  E- 3 

5.  "87  5  23  10E +3 

6.  1 96 1 04 "5 E  f3 

2.  3661567  1  E  +4 

46 


Table  VI- 1.  Nodal  Conditions  for  "Real  World"  Ephemeris  (Continued) 


Node 

Nodal  Crossing  Time 
/Mo/  Day/ Yr  \ 
\Hr/Min/Sec/ 

Position  XYZ 
(ft) 

Velocity  XYZ 
(fps) 

10 

3/ 1/80 

13/25/14. 56680 

-  1. 63910809E+7 
1.463  29058E+7 
-3. 17691 850E-4 

5.  96?  15899E+3 

6.  21640000E+3 

2.  36645896E+4 

1 1 

3/1/80 

14/54/41. 66063 

-  1. 64355506E+7 

1 . 45773503E  +  7 
-2. 74429 1 65E-4 

5.  94617  177E+3 

6.  23661260E+3 

2.  36677299E+4 

12 

3/1/80 

16/24/8.  03161 

-1. 6480367 6E+7 

1. 45214059E+7 
-  1. 76019997E-3 

5.  92479582E+3 

6.  25703989E+3 

2.  367027 24E+4 

13 

3/1/80 

17/53/33. 60059 

-  1. 65248092E+7 
1.44654088E+7 

1 . 48863  T 1 8E-4 

5.  90404820E+3 

6.  27813 694 E +3 

2.  367  27759E+4 

14 

3/1/80 

19/22/58.32904 

-  1. 65687151E+7 

1. 44090721E+7 

2.  9533  1619E-6 

5.  884  1  5436E+3 

6.  2992247 8E +3 

2.  36757447  E+4 

15 

3/1/80 

20/52/22. 17168 

-  1.661  227  50E  +7 

1 . 4  352590  IE  +7 
-2.  527 07 804 E- 4 

5.  8643  34  10E +3 

6.  3  2034 803 E +3 

2.  3 67 87 53  8E  t4 

16 

3/1/80 

22/21/45. 07135 

-  1 . 6655 1 805 E +7 

1 .  <*2959754  E  +-7 
-9.  013  14653K-4 

5.  844  3  2969E  +  3 

6.  34 1 30c  57  E  +3 

2.  36821387E  +4 

17 

3/1/80 

23/51/7.  1  1635 

-1. 669823  14E+7 
1.423961  18E+7 
-  1. 36532782E-3 

5 . 82390063 E +  3 

6 . 36 1 5  2469  E  +  3 

2 .  3  6  8  M  1  07  V.  +4 

18 

3/2/80 

1/20/28. 39917 

-  1 . 674053 88E  +7 
1.41829899E+7 

-  1 . ? 3762579E- 3 

8. 80307  0  17 E  t3 

6.  38107658E+3 

2.  3  0  889  3  68 E  +4 

Id 

3/2/80 

2/49/48.79290 

-  1. 67 827  64 9E +7 

1.41  258533  E+7 
-9. 0290061  IE-4 

5. 78157663  8:  *  3 

6.  40182983  8:  +  ) 

2.  3  6928  1408;  +4 

20 

3/2/80 

4/19/8. 28148 

-  1. 68255872E+7 

1 . 40682  145E+7 

-  9.  47885  1  95E-4 

5 .  7  5967  56  28:  +  3 
6.4  22077  19  8: +3 

2.  36904  6988:  +4 

47 


SECTION  VII 


COMPARISONS  OF  "MODEL"  AND  "REAL  WORLD"  EPHE  ME  RIDES 


A  trajectory  difference  run,  similar  to  those  described  in  Section  IV, 

was  made  for  the  model  and  the  real  world  ephemerides .  Figures  VH-la, 

b,  and  c  depict  the  radial,  in-track ,  and  cross-track  differences,  respectively. 

The  radial  difference  plot  shows  a  1-rev  periodicity  whose  amplitude  grows 

with  time  to  approximately  3300  ft  by  rev  20.  The  in-track  difference  plot 

also  shows  a  1-rev  periodicity  with  maximum  amplitude  about  5000  ft  and  a 

secular  term  which  grows  about  750  ft/rev.  Much  of  this  secular  growth 

can  be  attributed  to  the  differences  in  GM  between  the  two  geopotentials 

(the  model  GM  is  1 . 4076538841E+16  ft^/sec^,  while  the  real  world  GM  is 

3  2- 

1 . 4076468597E  +  1  6  ft  /sec  ).  The  cross-track  difference  plot  shows  a  clear 
1-rev  periodicity  and  modulation  at  some  multiple  number  ol  revs.  The 
maximum  amplitude  is  about  800  ft;  no  secular  growth  is  discernible. 

To  further  illustrate  the  differences  between  the  model  and  real  world 
ephemerides,  the  latitude,  longitude,  and  altitude  at  selected  perigees  through¬ 
out  the  20- rev  span  are  shown  in  Table  VII- 1  . 

From  the  foregoing,  it  is  reasonable  to  conclude  that  the  real  world 
ephemeris  is  in  a  sense  representative  of  the  real  world.  It  is  probably  noL 
as  complex  as  the  actual  real  world,  but  it  is  sufficiently  more  complex  than 
an  ephemeris  generated  by  conventional  models  for  the  purpose  .u  hand.  The 
difference  run  just  discussed  indicates  that  the  real  world  ephemeris  could 
not  readily  be  reproduced  by  conventional  modeling  techniques. 

The  content  and  format  of  the  final  data  tape  delivered  to  contractors 
is  discussed  in  Appendix  C. 
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Figure  Vll-ib.  Intrack  Difference  Between  "Model"  Fphemeris 
and  "Real  World"  Fphemeris 
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Figure  VH-ic.  Crosstrack  Difference  Between  "Model" 
Fphemeris  and  "Real  World"  Fphemeris 


Figure  Vll-ic.  Crosstrack  Difference  Between  "Model" 
Ephemeris  rind  "Re a  1  World"  Ephomeris 


Table  VII- 1.  Selected  Perigees  from  "Model'1  and  "Real  World" 
Epheme  rides 


Perigee 

i  '  "  -  11 

"Model" 

"Real  World" 

1 

LAT 

45.  17868651° 

LAT 

45. 16178965° 

LONG 

193. 33637181° 

LONG 

193. 32125882° 

ALT 

81.36217  n  mi 

ALT 

81.3  6947  n  mi 

5 

LAT 

45.67476127° 

LAT 

45.72686973° 

LONG 

104.69701390° 

LONG 

104.74  014037  ° 

ALT 

81 . 36473  n  mi 

ALT 

81.4619°  n  mi 

10 

LAT 

46. 38336256° 

LAT 

46.44256510° 

LONG 

354. 06792672° 

LONG 

354. 1 1780809° 

ALT 

81. 16635  n  mi 

ALT 

j 

81.44194  n  mi 

15 

LAT 

46. 83655544° 

LAI 

4  6.  7  8  5  7  ~  5 1  >  'i  ° 

LONG 

243. 30253604° 

i.o.v ; 

24  !.  254  0  54  21  ° 

ALT 

80 . 7  084 1  n  mi 

A 1 .  i 

81.  1  6  Ri.l  n  mi 

20 

LAT 

47. 4 84 27 "83° 

l.A  1 

47 . 4 '< 200754° 

LONG 

13  2.  816907  27° 

I  .<  >NG 

1  52.  7  0  3  7  0  <  1 0  *  1 0 

ALT 

80.  54  5  62  n  mi 

A  Li 

81.  1  24  60  n  mi 

LATITUDE 


LONGITUDE 


t  iiMirt*  A-l.  Ground  Trace,  Revs  !  through  5 


LONGITUDE 


6 


Figure  A -2.  Ground  Trace,  Revs  f>  through  10 
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APPENDIX  C 


'REAL.  WORLD  '  E  PH  EM  ERR  TAPI- 

and  its  format 


The  real  world  data  tape  was  generated  on  Control  Data  Corporation 
(CDC)  6000  series  machines.  It  is  an  unlabeled  binary  tape  written  in  odd 
parity  at  a  density  of  556  bpi. 

A  seven-track  non- return-to-zero  (change-on-ones )  recording  scheme 
is  used.  Magnetic  particles  on  the  tape  are  aligned  in  either  the  positive  or 
negative  direction.  A  binary  "l'1  is  recorded  by  reversing  the  alignment 
(polarity);  no  polarity  reversal  results  in  a  "0.” 

A  frame  of  tape  data  consists  of  one  6-bit  data  character  and  one  parity 
(check)  bit  for  each  character.  Tracks  0  through  5  specify  the  characters 
while  track  6  holds  the  parity  bits  (see  Figure  C-l).  In  binary 


800  BPI  =0  0023 
336  BPI  =  0  0016(00036) 
200  BPI  =0  005(0.01) 


Pr 

J* —  I  ,  ... 


r 


A 


(READ)  0032 


TRACK 

DETAIL 


0  048 


-REC0R0- 


i- 


NOTE 


OXIDE  SIOE  UP  ON  DIAGRAM,  RECORDING 
HEAD  ON  SAME  SIDE  AS  0XI0E 


200  BPI  =  0  02(0  04) 

556  BPI  =  0  0074(0  0148) 
800 BPI =  001 


Figure  C-l.  Hit  Assigiments  on  Tape 

format,  the  parity  bit  is  chosen  so  that  the  total  number  ot  1  '  bits  in  a  in¬ 
line  is  odd.  Additional  information  on  I  lie  CDC.  tape  processing  hardware 
and  procedures  is  available  in  CDC  Publication  \’o.  60156  100. 
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CDC  employs  a  60- bit  floating  point  word;  its  format  is  shown  in 
Figure  C-Z.  Floating  point  arithmetic  takes  advantage  of  the  ability  to 

COEFFICIENT  6IASE0  INTEGER 

sign  exponent  coefficient 

H~  n  j  <3~~]  i 

59  58  46  47 

Figure  C-2.  Floating  Point  Format 

express  a  number  as  kBn,  where  k  -  coefficient,  B  =  base  number,  and 
n  -  exponent.  The  base  number  is  constant  (2)  for  binary-coded  quantities 
and  is  not  included  in  the  general  format.  The  binary  point  is  considered  to 
be  to  the  right  of  the  coefficient,  thereby  providing  a  48-bit  integer  coefficient, 
the  equivalent  of  about  14  decimal  digits.  The  sign  of  the  coefficient  is  car¬ 
ried  in  the  highest  order  bit  of  the  packed  word.  Negative  numbers  are 
expressed  in  one's  complement  notation.  The  11 -bit  exponent  carries  a  bias 
of  21(^(2000g)  when  packed  in  the  floating  point  word.  The  CDC  procedure  used 
in  floating  point  arithmetic  is  described  in  CDC  Publication  No.  60100000. 

The  tape  contains  one  file  which  consists  of  216  records.  Each  record 
contains  501  words;  each  word  is  60  bits. 

The  first  word  in  each  record  is  a  bookkeeping  word;  it  is  an  integer 
quantity  equal  to  the  number  of  data  words  in  that  record.  Since  there  are  no 
short  records  on  this  tape,  this  word  is  always  the  integer  500.  The  remain¬ 
ing  500  words  are  segmented  into  50  data  frames  of  1  0  words  each.  Table  C-l 
gives  the  ordering  and  units  of  the  10  data  words. 

All  position,  velocity,  and  acceleration  components  are  given  in  an 
i- ;  ( ;  J  coordinate  frame.  The  X  axis  of  this  coordinate  system  is  pointed 
toward  the  mean  equinox  of  date  at  0  hr  GMT  of  the  date  of  epoch.  The  X 
and  V  axes  form  the  plane  of  the  true  equator  at  epoch.  The  V.  axis  is 
perpendicular  to  the  true  equator  at  epoch.  'The  set  is  right  handed,  with 
V  lying  east  of  X  and  /  being  positive  north. 


0 

BINARY 

POINT 
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Table  C-l.  Data  Frame  Format 


Word 

Quantity 

Description  and  Units 

1 

time 

Time,  sec  from  epoch 

2 

X 

3 

Y 

ECI  vehicle  position,  it 

4 

Z 

5 

X 

6 

Y 

ECI  vehicle  velocity,  fps 

7 

Z 

8 

x_ 

D 

9 

VD 

2 

k  ECI  total  atmospheric  drag  acceleration,  ft/sec 

10 

•  • 

Z  D 

The  data  rate  on  the  tape  is  one  data  set  for  each  10  sec.  The  first  time 
on  the  tape  is  0  sec;  it  corresponds  to  the  epoch  time  of  22  hr  30  min  0  sec 
GMT  on  29  February  1980.  The  last  time  on  the  tape  is  107990  sec. 


C-3 


REFERENCES 


1 .  C .  C .  Tonies,  et  al.  ,  TRACE  Orbit  Determination  Program  -  Version  D, 
TR-669(9990)- 3 ,  The  Aerospace  Corporation  (September  1  q66) . 

2.  G.  Buechler  and  D.  Walker,  TRACE  66  Orbit  Determination  Program  - 
Volume  III;  Trajectory  Generation  Equations  and  Methods, 

TOR -0066(9320) -2,  Vol.  Ill,  The  Aerospace  Corporation  (25  April  1970). 

3.  W.  H.  Guier  and  R.  R.  Newton,  "The  Earth's  Gravity  Field  as  Deduced 
from  the  Doppler  Tracking  of  Five  Satellites,  "  J.  Geophys.  Res.  ,  70, 
4613-4626  (September  1965). 

4.  R.  W.  Bruce,  An  Atmosphere  Density  Model  Recommended  for  Analysis 
of  Low  Altitude  Satellite  Orbits,  TOR- 1 001  (2 1 1  0-0 1 ) -8,  The  Aerospace 
Corporation  (30  August  1966). 

5.  W.  M.  Kaula,  "Determination  of  the  Earth's  Gravitational  Field," 

Rev,  of  Geoph. ,  J_  (4)  (November  1963). 


R  - 1 


